Switchable Quantum Anomalous Hall state in a strongly frustrated lattice magnet 
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We establish that the interplay of itinerant fermions with localized magnetic moments on a checkerboard 
lattice leads to magnetic flux-phases. For weak itineracy the flux-phase is coplanar and the electronic dispersion 
takes the shape of graphene-like Dirac fermions. Stronger itineracy drives the formation of a non-coplanar, 
chiral flux-phase, in which the Dirac fermions acquire a topological mass that is proportional to a ferromagnetic 
spin polarization. Consequently the system self-organizes into a ferromagnetic Quantum Anomalous Hall state 
in which the direction of its dissipationless edge-currents can be switched by an applied magnetic field. 
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Introduction. — The study of topologically non-trivial states 
of matter is one of the hottest topics in present day condensed 
matter physics. An understanding of topological states re- 
quires a theoretical paradigm that goes far beyond the concept 
of global symmetry breaking that has originally been laid out 
by Landau. It is remarkable that the theoretical predictions 
on the existence of various topologically ordered states have 
rather swiftly led to the discovery of an entirely new class of 
materials, the topological insulators IJ\^. Recent pioneer- 
ing experiments have confirmed the key signatures of non- 
trivial topology in certain materials, e.g. spin-momentum- 
locked undoubled Dirac fermions |5 -7 1 and the Quantum Spin 
Hall (QSH) effect |8|. These topological insulators are time- 
reversal (TR) invariant generalizations of the first, much older, 
topological state of matter, the famous Integer Quantum Hall 
states |9 10| that are induced by a magnetic field, which ob- 
viously breaks TR symmetry. 

In a seminal work in 1988, Haldane established that a mag- 
netic field is not required to induce states with the same topol- 
ogy as IQH states ifTTIl . It was shown that adding complex 
hopping to a graphene-like Hamiltonian for electrons on a 
honeycomb lattice opens up topologically nontrivial gaps at 
the Dirac points, which yields a topologically ordered, insu- 
lating state, refeiTed to as a Quantum Anomalous Hall (QAH) 
state. An important feature of QAH states are edge channels, 
in which cuiTent can run only in one direction; in contrast to 
QSH states, on a single edge the opposite spin channel carry- 
ing the opposite current is absent 1 12 1. QAH states would thus 
allow very robust, dissipationless charge transport along edge 
channels, as backscattering would be completely suppressed. 
However, while signatures of QAH behavior have been re- 
ported in some compounds jTsUTSi . the QAH state is the only 
one among these topologically insulating states that remains 
to be unambiguously identified in experiment. 

The experimental difficulty is mirrored by the frailty of the- 
oretical mass-generating mechanisms for a graphene-like ki- 
netic energy with a linear dispersion at the Fermi level. TR- 
symmetry breaking via (magnetic) order requires rather spe- 
cific and strong longer-range Coulomb interactions |16|, be- 
cause the Dirac cones' vanishing density of states at the Fermi 
level renders interaction-driven ordered states energetically 
less favorable. QAH states can more readily be induced in 



models with a finite density of states lfT7l - [T9]| . especially in 
cases of quadratic band crossings ||20ll , as for instance found 
in the checkerboard lattice, which exhibit a weak-coupling 
instability |20-22|. Another approach has been to consider 
spin-orbit coupled magnetic semi-conductors ll23l or spin- 
polarized QSH states ifTH . 

Our starting point is instead the Kondo lattice model, which 
provides the most general context for the study of the interplay 
between localized spins and itinerant electrons. We will show 
that, depending on parameters, its ground state on the checker- 
board lattice can feature massless Dirac cones or a chiral QAH 
state. Moreover, because the spin texture underlying the QAH 
state has a net ferromagnetic (FM) moment, we have a di- 
rect switching mechanism between ground-states with differ- 
ent chirality just by flipping the FM polarization. The ensuing 
possibility to reverse the direction of an edge current by an 
external magnetic field is an attractive feature in the context 
of spintronics. 

In particular we focus on the case of electrons strongly cou- 
pled to large localized spins that in turn interact via a strongly 
frustrated antiferromagnetic superexchange. In the absence 
of charge carriers, the magnetic interactions give rise to a 
highly degenerate ground state manifold (GSM) comprised 
of all spin configurations that obey certain local constraints, 
which go by the name of spin-ice rules, as a reference to three- 
dimensional spin systems that remain disordered down to the 
lowest temperatures like the H protons in water ice [24|. It is 
clear that in a frustrated spin-systems governed by ice rules, 
doping of itinerant charge carriers will (partly) lift the macro- 
scopic the degeneracy of the GSM |25 1 because of the kinetic 
energy competing with the local ice-rule constraints. We will 
show that for not too large kinetic energies, a unique flux- 
phase ground-state is selected in our case, whose low-energy 
states are described by a massless Dirac equation. When the 
kinetic energy becomes stronger, enhancing the competition, 
these Dirac fermions acquire a topological mass, and the fer- 
romagnetic QAH state emerges. 

Model and method. — We thus consider itinerant electrons 
coupled to localized core spins, described by the canonical 
one-band double-exchange model with a competing antiferro- 
magnetic (AFM) super-exchange interaction, on the checker- 
board lattice at the density of one electron per two sites, where 



2 



•A.B jVx 
















-(p 






-(p 




-<f 






-cp 





(c) 



FIG. 1. (Color online) (a) Shows the localized spin texture of 
the "flux" phase; all spins are in plane, for number labeling see 
text, (b) Representation of the spins in the "umbrella" spin-chiral 
state, as obtained with MCMC+optimization for a 16x16 lattice with 
Jaf = 0.105, where S = 0.148 compared to 0.141 as would be ex- 
pected analytically; all spins have been translated to a single site, (c) 
Staggered flux arrangement on the checkerboard lattice correspond- 
ing to this insulating umbrella phase. Red arrows indicate the gauge 
choice of the "gauge flux". Note that the indicated flux threads the 
squares, so half of it threads the triangles. 



the kinetic energy promoting ferromagnetic (FM) spin corre- 
lations is strongest. The resulting Hamiltonian is 



H = -Y1 ^^MUj + H.c.) + Jaf ^ S, 



(1) 



where ipl (ip- ) creates (annihilates) a fermion on site i. Here 
we have assumed for simplicity that an infinite Hund's rule 
perfectly aligns the fermion spin to the localized spins S^, 
but we have verified that the results presented in this Let- 
ter remain valid also for large but finite coupling. As the 
onsite spins are classical and one can take without loss of 
generality |Si| — 1, they are completely specified by polar 
and azimuthal angles {9i, 4^). As hopping and superexchange 
are identical along the "straight" and "diagonal" edges of the 
checkerboard lattice, see Fig. [TJa), both types of bonds are 
included in the sum over (ij). The hopping amplitude de- 
pends on the core spins as tij = t[cos{6i/2) cos{9j /2) 



i(0,/2) sm{9j/2) 



')] ||26|. AFM super-exchange is 



given by Jaf and all energies will be measured in units of 
the hopping amplitude t. We use Markov Chain Monte Carlo 
(MCMC) simulations to treat the classical spins, where the 
weight of a spin configuration is given by the free energy of 
the effective fermionic Hamiltonian, as obtained by exact di- 
agonalization |f26l. We have performed calculations on lat- 



tices with N = 



12^, 16^, and 20^^ sites. MCMC calcula- 



tions were supplemented with an energy optimization in order 
to suppress thermal fluctuations 1.27 J . 

Checkerboard lattice magnet. — Let us first discuss the two 
terms in the Hamiltonian separately, coiTesponding to the lim- 
its Jaf and Jaf oo, and turn to them compet- 
ing at the next stage. For vanishing super-exchange cou- 
pling Jaf 0, one finds a saturated FM state, equivalent 
to free spinless fermions on a checkerboard lattice. There are 
two bands 
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2t and E_ = —2t — At cos cos ky and 
the density of states (DOS), D{uj) = {j^ J2k ~ '^fe))' 
shown in Fig. |2|a). The first Brillouin zone (BZ) is given by 
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FIG. 2. (Color online) Electronic density of states of (a) the FM 
(spinless) phase, (b) the "flux" phase (in red) and spin-chiral um- 
brella phase for various S (see text). 



y) G {\kx + ky\ < tt} n {\kx — ky \ < vr}. The Heisen- 
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berg term dominating for Jaf — oo is geometrically frus- 
trated and supports an infinitely large classical degeneracy of 
spin ground states. This can be seen by noting that it can be 
rewritten in terms of the total spins on the crossed plaquettes 
Sp = Si + S2 + S3 + S4, giving 



Jaf ^ Si 



Jaf 



^ Sp • Sp - JafN, (2) 



where N is the number of sites. Clearly, the lowest energy is 
obtained when S-p = 0. This local requirement on the spins 
of each crossed square is similar to the spin-ice rule in the 
pyrochlore lattice. For Ising spins |28| the rule coiTesponds 
to "two up-two down". If noncollinear spin arrangements are 
permitted, as we consider here, the class of states S-p = is 
further increased. 

For dominant super-exchange coupling Jaf ^ 1, the mag- 
netic order is expected to belong to the highly degenerate 
ground-state manifold (GSM) fulfilling S-p = 0. The ki- 
netic energy can leave the degeneracy intact, lift it partially or 
lift it completely, singling out a unique non-degenerate spin 
arrangement. Our MCMC calculations show that the latter 
is the case, the electrons pick out a particular coplanar, but 
not collinear, state that is schematically depicted in Fig. [TJ^)- 
Non-diagonal bonds connect orthogonal spins, while diagonal 
bonds connect AFM spins, effectively excluding them from 
the hopping term. Going around a square plaquette, the elec- 
trons pick up a phase e", corresponding to a time-reversal in- 
variant flux of TT, and this special "flux" phase has been shown 
to arise in models for high-Tc superconductors Il29ll30l and in 
the double-exchange models on the square lattice ||3Tti34 1. On 
the unfrustrated square lattice, it competes with the Neel state 
for strong Jaf li31ri33J . but since it fulfills the ice-rules, it 
remains stable for Jaf > 0.12 on the checkerboard lattice. 

The DOS of the flux phase shows semi-metallic behaviour 
[see Fig. |2|b)] that originates from two Dirac points in the 
spectrum. Low-energy excitations are described by a rela- 
tivistic Dirac equation, in full analogy with graphene |35|. 
The core-spin texture = (0^,0^) of the flux phase can 
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be written as = (7r/2, (i - l)7r/2) with i = 1,2,3,4 
[see Fig. [TJa)]. Even though the magnetic texture has a 4- 
site unit cell, the two-site electronic unit cell need not be en- 
larged. The electronic Hamiltonian, in the {ipj^jip^g) basis, 
is then given by i?(k) = d(k) • cr, with he Pauli matrices 
cr = (fj^, cr^, fj^) and d(k) = — (cos fcj; + cos A:^, cos fc^: — 
cosfcyjO). The band structure of this state is shown in 
Fig. [3ja). The two inequivalent Dirac points, or valleys, are 
located at Mj- = (±7r/2, 7r/2). Expanding around the Dirac 
points yields an effective low-energy Hamiltonian with the 
two Dirac spinors = {^p\{M+ + p),'iJj^g{M+ + p))and 

*2 = (V'V^^- + P);V'L(^- + P))- Rotations on mo- 
mentum and spinor components bring us to the familiar form 
H{q) — ® {Qx^'^ + ^), where i/-Pauli matrices act 
on the valley index. This is equivalent to graphene, with two 
valleys around which the electrons are described by the Dirac 
equation and an interesting mapping exists between graphene 
and the flux phase |f36l . One important difference to graphene 
is that we have here no spin degeneracy, as the spin degree of 
freedom was integrated out by tying the fermion spin to the 
core spins. 

Massive QAH Dime fermions. — Having established that a 
not-too-large electronic kinetic energy selects a unique non- 
collinear pattern for the checkerboard double-exchange mag- 
net, which has a graphene-like Dirac spectrum, we consider 
next what happens upon an increase of the itineracy. Low- 
ering Jaf, we find that the magnetic interactions enforc- 
ing the tetrahedron rules are overcome by the electronic ki- 
netic energy for Jaf ^ 0.12. The transition is continuous 
and can be understood as a tilting of the flux-pattern out of 
the plane, forming an "umbrella". An example is shown in 
Fig. [TJb): the spins fall along four directions, whose pro- 
jections onto the x-y plane mirror the "flux"-phase pattern, 
but there is an additional FM component along the z axis. 
The spins can be described using an Ising variable s = ±1 
(which will turn out to correspond to a scalar spin chiral- 
ity) and a continuous parameter 5 giving the tilting along ^z: 
{Af (5)} = (7r/2 + 5, s{i - l)7r/2), where i = 1, 2, 3, 4 again 
runs around a crossed plaquette. A similar scenario, but in- 
volving a more complex 8-site unit cell and leading to Chern 
numbers ±2, arises on a square lattice with longer-range cou- 
plings when nearest-neighbor hoppings are strongly modu- 
lated i34;|. 

The scalar spin chirality of the state is defined as x — 
'YliT Si ■ ^ Sfc' where the sum is over all triangles T of the 
checkerboard lattice, and • Sj x is taken in the counter- 
clockwise direction. The chirality as function of 5 is plotted 
in the inset of Figj2|b) for umbrella states A*, it is x ~ ~s5 
for small 5. The label ±s decides the sign of the chirality for 
(5 > and is related to a (counter-)clockwise rotation of the 
spin projection onto the x-y plane. The umbrella states, in ad- 
dition to a continuous spin rotation symmetry, thus also break 
a discrete Z2 symmetry. Since breaking a discrete symmetry 
in 2D is possible at finite temperature, this opens up the pos- 
sibility for chiral ordering in absence of long-range magnetic 
ordering li37J . 



The effect of the tilting on the electronic degrees of free- 
dom is to break time-reversal symmetry, as fluxes through ele- 
mentary plaquettes are related to the solid angle subtended by 
the spins surrounding the plaquette. Calculating the hoppings 
in the umbrella states, we find that hopping on the straight 
bonds is given by t\ = e^**'^/^(l — sisu\6) / ^/2, with \t\\ = 

y^(l + sin^ 5)/2 = ti and 0^ ~ arctan(— s sin 5) — s7r/4 = 
[see Fig. [TJc)]. In addition, hopping along the diagonal 
bonds is no longer but t2 = — sin (5, independent of chiral- 
ity. Using these expressions we write the effective Hamilto- 
nian for the electrons as 

JJ'*(k) = d(k) • cr + d°(k)cr°, with (3) 
(i°(k) — ~2t2 cos kx cos ky, (i'^(k) = 2^2 sin fc^; sin fcj^ 
(i^(k) — —2ti cos (j)" {cos kx + cos ky) , and 
(P(k) — — 2ti sin0''(— cos kx + cos ky), 

where the two states referred to by the matrices are again the 
two sites of the unit cell and (7° is the unit matrix so that the 
Dirac Hamiltonian above is recovered for 6 — 0, implying 
(f)" =7r/A and = 0. From the DOS [Fig.gb)] and the band 
structure [Fig. [3]; a)], it is clear that finite S ^ opens a gap for 
the Dirac cones. Since the hoppings are complex and the diag- 
onal bonds have been activated, both time-reversal and parity 
symmetries are broken, allowing a QAH state ll38l . To estab- 
lish that the gapped state is indeed topologically non-trivial, 
we calculate the Chern number C„ = §93 z ' -^C^)' 
where A(k) = (nk|Vk|?^k) is the Berry connection and find 
C = sgn(t2)sgn(sin 20''). Chirality and Chern number hence 
perfectly correlate and we observe that inverting the magnetic 
polarization S ^ —S flips both the spin chirality and the 
Chern number. The off-diagonal Hall conductivity as a func- 
tion of chemical potential, obtained from Eq. ([3]) for S = 0.3, 
is shown in Fig. |3jb). Figure [3|c-e) shows the effect of non- 
trivial topology on the edge of the system: chiral edge states 
connect valence and conduction band. As can be seen by com- 
paring Figs.|3jd) and[3|e), the direction of the edge currents 
can be reversed by inverting the spin chirality. The latter can 
be easily manipulated by a small magnetic field. 

The observation that spin configurations of the umbrella 
states are continuously connected to the coplanar flux phase 
suggests that the electronic QAH state can be understood from 
the low-energy physics at the Dirac points. We will demon- 
strate this by analyzing, in the spirit of Ref [11], the system 
in a presence of an external magnetic field B and then take 
the limit of S — > 0. Focusing on the low-energy theory, Sk — 
k — M-y (7 — ±), we introduce the magnetic field by way of a 
Peierls substitution hSk it, where tt is the dynamical mo- 
mentum whose components satisfy the commutation relation 
[rrxjTry] = leBh. We obtain two independent Hamiltonians 
for the two Dirac points, H-y = vpiTt^yO-^' + '^^cr^) + m^a^, 
which indeed has the appearance of the relativistic Dirac equa- 
tion in a magnetic field. As can be seen by comparing to 
Eq. Q, our mass term = 2'^t2 is a direct consequence 
of finite t2 — ~ sin 5, and thus of finite chirality 5 ^ {). Op- 
erators TT^ and TT^ are derived from tt and satisfy the com- 
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FIG. 3. (Color online) (a) The band structure of the gapless flux 
phase (red) and the insulating chiral phase (blue, S — 0.3) along a 
path in the Brillouin zone specified in the inset, (b) Quantized Hall 
conductivity in the chiral state {5 = 0.3), when the Fermi level is 
in the gap, the quantized value depends on the chirality of the spin 
state. The inset shows the calculated chirality of the states A*, where 
dashed (solid) corresponds to + (— ). (c-e) Spectrum of the flux phase 
calculated for a strip geometry, which explicitly shows the edge states 
at the open boundary, (c) yr-flux phase (5 — 0.0) exhibits edge states 
similar to graphene. (d,e) Chiral gapped phase (5 = 0.2); chiral edge 
states connect valence and conduction bands. The states drawn with 
solid (dashed) lines lives on the top (bottom) edge. The chirality in 
(e) is reversed with respect to (d), the right- and left-moving states 
are consequently exchanged. 

mutation relation [7riJ,,7r^] = —ijsm{2(j>'')eBh. Relativis- 
tic Dirac fermions in a magnetic field are known to exhibit 
zero modes in their spectrum [39J, which cause the charge 
density imbalance in the ground state, potentially leading to 
an integer QAH effect. Here, the zero modes have energy 
_Bo,7 = — 7m.7Sgn(sin 2(/>*)sgn(ei?), and the spectrum is 
asymmetric when m+ and m_ have opposite sign. Following 
Haldane, we obtain the off-diagonal conductivity in the limit 
B — 0, (Txy — ve^/h, where — ^sgn(sin2(/)*)[sgn(m+) — 
sgn(m_)] = sgn(i2)sgn(sin20*). Hence, the gapped QAH 
umbrella state can be interpreted as Dirac fermions becom- 
ing massive, with masses of opposite sign, indeed, the (k) 
component of Eq. ([3]) has opposite sign at the two Dirac points 
(±7r/2, 7r/2). A sublattice potential, which also gaps out the 
Dirac fermions, would in contrast lead to equal masses, and 
the edge states would not cross the chemical potential. 

Discussion and conclusions. — We investigated the inter- 
play of itinerant electrons with a frustrated AFM spin back- 
ground on the checkerboard lattice using Monte-Carlo meth- 
ods. From the macroscopically degenerate spin-ice ground- 
state manifold, which optimizes the AFM interactions, the 
electron kinetic energy selects a unique magnetic ground state. 
The low-energy electronic states of the selected tt-Aux phase 
are given by a relativistic Dirac equation, and slightly stronger 



kinetic energy induces a spin chirality, from which the Dirac 
fermions inherit a topologically nontrivial mass. The Kondo- 
lattice model on the checkerboard model thus provides a di- 
rect realization of Haldane's proposal for obtaining a QAH 
state 1 11 1. In addition, the QAH state's chirality is coupled to 
a FM spin polarization and the direction of the edge currents 
can thus be switched by a magnetic field, an alluring property 
for quantum spintronics applications. 

Our findings are also relevant in the broader context of 
fractionalization of quantum numbers in two different ways. 
First, the QAH state on the checkerboard lattice is actually a 
prominent candidates for hosting a fractional quantum-Hall- 
like state without a magnetic field, because the topologically 
nontrivial band can be made almost flat by tuning hoppings 
and flux |21 , 40, 41). It turns out that even if the flux and 
effective hoppings ti and ^2 that emerge in the "umbrella" 
configuration do not lead to very flat bands, additional longer- 
range hopping — 2t3(cos 2/0^; + cos2fcj,) can give a ratio of 
band gap vs. band width of w 5 for 5 — 0.3, considerably 
less than ratios achievable by tuning all parameters 121] or in 
^2g-orbital systems [42, 43 1, but comparable to Cg [42] sys- 
tems or a square-lattice model ll44l . Second, it was recently 
demonstrated that vortex defects of the localized magnetic or- 
der underlying a QAH state can carry fractional charge and 
spin quantum numbers in the electronic sector P3]| . 
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